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This paper discusses the characteristics, and some of the statistics of the general climatological 
aspects of considering atmospheric oxygen and water vapor absorption (the two principal contributors 
to atmospheric gaseous absorption in the 6 to 45 Gc/s frequency range) to be decreasing exponential 
functions of height from their surface value. The paper is intended primarily as an examination of 
these climatological aspects and the maps included are not primarily intended for application. The 
dry (oxygen) term of this biexponential model is extremely well behaved, but the wet (water vapor) 
term depends on time of year, geographic location, air mass, and general climatology of the areas for 
which application is intended. 

Gaseous absorption of radio energy in the 6 to 45 Gc/s frequency range arises principally from 
water vapor and oxygen. The combined absorption may be represented by a biexponential decrease 
from the surface values. The oxygen absorption is well described by the model, as is water vapor 
absorption except in regions of very low water vapor content. 



1. Introduction 

In the frequency range 6 to 45 Gc/s, the atmospheric, 
gaseous absorption, y, at a particular frequency, ex- 
pressed in dB/km, arises primarily from oxygen (dry) 
absorption, yd, and water vapor (wet) absorption, y w ; 
i.e., 



y = ya + yw 



(1) 



Zhevankin and Troitskii [1959] have indicated that 
ya and y w can be represented as exponential functions 
of height, z, above the earth's surface, 



yd = y<i,o exp I— j, y w = yw,o exp ( — ), 



(2) 



where ya, o and y w , o are the values of ya and y w , re- 
spectively, at the earth's surface, and Hd and H w are 
called the "scale heights" of yd and y w - This model 
is known as the "biexponential" model of absorption. 
The values of yd and y w are often called the "dry" 
and "wet" terms of y. 

This paper primarily examines the possible extent 
to which this model can be employed on a worldwide 



1 This paper was presented at the World Conference on Radio Meteorology, Sept. 14-18, 
1964, National Bureau of Standards, Boulder, Colo. 



climatological basis, to either support or not support 
application of the model at a given location. Since 
the data used in this paper are not extensive, the paper 
is not primarily intended to be used in applications. 
More extensive data taken during the seasons dis- 
cussed are advisable to meet this end. 

The scale height for the dry term in the 6 to 45 Gc/s 
frequency range can be written as 



H d = 



To 



c + ba 



(3) 



where T is the surface absolute temperature, a is the 
temperature lapse rate with height, which is assumed 
to be a constant for the entire atmosphere, and b and 
c are constants determined from thermodynamic 
considerations. 

2. Statistical Fitting of the Biexponential 
Model 

Radiosonde data for 49 climatically and geographic- 
ally diverse stations were fit, by means of a least 
squares regression method, to the biexponential 
absorption model. Two radiosonde profiles for each 
station were chosen, one in the month of February, 
the other in the month of August. Each profile was 
so chosen as to be reasonably representative of average 
wintertime and summertime meteorological conditions 
at the site where the profile was obtained. 
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Figure 1. 



Exponential fits to oxygen and water vapor absorption 
versus height at Balboa, Canal Zone. 



Table 2. Fit of y w = y w , exp [- z/H w ] 



Frequency 



Surface water 
vapor absorption 



Scale height 



Correlation 

coefficient of 

fit 



Standard error 

of estimate of 

y w , S. E. 



at Barrow, Alaska, February 20, 1952, 0500 UT 



Gc/s 
6.0 

22.406 
42.0 


dB/km 
0.00006 
.00498 
.00245 


km 

-2.7747 
-2.3254 
-2.7792 


+ 0.64063 
+ .73546 
+ .63995 


dB/km 
±0.00277 
± .62074 
±.11748 



at Sopporo, Japan, August 26, 1952, 1500 UT 



6.0 

22.406 
42.0 



0.00418 
.48160 
.17860 



2.4989 
4.1971 
2.4961 



-0.97903 

- .95934 

- .97912 



at Balboa, Canal Zone, August 16, 1952, 1500 UT 



6.0 

22.406 
42.0 




-0.99832 

- .99256 

- .99833 



at Denver, Colo., U.S.A., August 16, 1952, 0300 UT 



±0.00018 
± .03513 
± .00777 



±0.00007 
± .02475 
± .00300 




at Guryev, U.S.S.R., August 23, 1957, 0000 UT 



Table 1. Fit of y d = y (i ,o exp[— z/H d ] 



Frequency 



Surface oxygen 
absorption 



Scale height 



H d 



Correlation 

coefficient of 

fit 



Standard error 

of estimate of 

yd, S. E. 



at Barrow, Alaska, February 20, 1952, 0500 UT 



Gc/s 


dB/km 


km 




dB/km 


6.0 


0.01167 


3.5616 


-0.99473 


±0.00032 


22.406 


.02374 


3.5387 


- .99463 


±.00067 


42.0 


.17433 


3.5568 


- .99471 


±.00486 



at Sopporo, Japan, August 26, 1952, 1500 UT 




at Guryev, U.S.S.R., August 23, 1957, 0000 UT 



6.0 
22.406 
42.0 




-0.99161 
- .99152 
-.99159 



±0.00024 
±.00049 
± .00362 



6.0 

22.406 
42.0 



0.00262 
.29374 
.11190 



2.5987 
4.1792 
2.5952 



-0.74450 
- .65330 
-.74489 



±0.00033 
± .04856 
± .01420 



Figure 1 shows an example of one of these profiles 
for Balboa, Canal Zone, on a typical August day. On 
figure 1 are plotted the data points of both water vapor 
and oxygen absorption determined from the meteoro- 
logical parameters obtained during the radiosonde 
ascent. Also plotted on figure 1 are the exponential 
fits [as in (2)] resulting from least squares regression 
of the data in the figure. Tables 1 and 2 also show 
some of the data for exponential fits of oxygen and 
water vapor absorption. 2 Note that Barrow, Alaska, 
in table 2, a poor fit for the water vapor absorption, 
has negative H w 's and positive correlation coefficients. 
This implies that y w in this case is actually better 
represented by an increasing exponential function of 
height rather than by a decreasing exponential func- 
tion, and points out a phenomenon peculiar to arctic 
wintertime conditions. 

Figures 2 and 3 show Ha and H w as functions of fre- 
quency at a near arctic (Verkhoyansk, U.S.S.R.) loca- 
tion and at a tropical (Balboa, Canal Zone) location. 
The spectra of Ha and H w for the two stations plotted 
in figures 2 and 3 were obtained from two least squares 
regression fits. Each of the two regression fits were 
from radiosonde data composing a given single profile 
which was used as representative of seasonal condi- 



2 A forthcoming NBS Technical Note on the material covered in this paper will include 
the complete data of the 49 station sample analogous to tables 1 and 2. These tables will 
include regressions for many more frequencies in the 6-45 Gc/s range. 
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tions at the station. Each plot shows the variation 
of Ha and H w in the location from summer (August) 
to winter (February). It is important to note the con- 
stancy with respect to frequency of the Hd curves even 
at the water vapor resonance frequency, 22.4 Gc/s. 
Thus, a single value of Ha determined from (3) may be 
used to represent Ha in the 6—45 Gc/s range, but a 
different value of H w is required for each frequency. 

Figures 4, 5, 6, and 7 are plots of Ha versus ya, o and 
H w versus y w , o using as data the results of the least 
squares biexponential fits discussed above on the 49 
stations. These figures also include the wintertime 
(February) and the summertime (August) values on 
the same chart. Figure 4 shows Ha versus ya, o at a 
wavelength, X, of 3 cm — well off water vapor reson- 
ance, and figure 5 shows Ha versus ka, o at water vapor 
resonance, A. = 1.33 cm. 

A least squares fit to the data of figures 4 and 5 of 
the form: 



Ha = Bya,o + A,(A and B are constants) 



(4) 



gives reasonably high correlation coefficients, r, of 
— 0.78 and —0.76 as seen on the figures. This is 
understandable, on the basis of (3), which strongly 
indicates the dependence of Ha on surface conditions. 
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Figure 2. Dry and wet term scale heights versus frequency at 
Verkhoyansk, U.S.S.R. 
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However, the data do not support the strong possibility 
of the similar expression, 

H w = Cy w , o + £), (C and D are constants) (5) 

for water vapor, as seen on figures 6 and 7. 
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Figure 4. H d versus 7d, o away from the y = 1.33 cm water vapor 
resonance. 
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FIGURE 3. Dry and wet term scale heights versus frequency at FIGURE 5. H d versus y d . o near the y=1.33 cm water vapor 
Balboa, Canal Zone. resonance. 
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Table 3. Values of actual total integrated absorption, T a = yds, and exponentially computed values T exp = T d + T w where r d = 

Jo 

exp (—J dz, t w = I y w , exp f— J dz, d = 90° 



ya.o 





Balboa, Canal Zone (August) 


Sopporo, Japan (August) 




Barrow, 


\Jaska (February) 






Ta 


T exp 


To 


T„. 


8t=t„— T exp 


Ta 


Texp 


Trf 


T u - 


8t 


T« 


Texp 


Trf 


Tw 


8r 


Gels 
6.0 
22.406 
42.0 


dB 

0. 043 
2.119 
0.949 


dB 

0.045 
2.321 
0.988 


dB 

0.034 
.069 
.511 


dB 

0.011 

2.252 

0.477 


dB 
-0.002 
-.202 
-.039 


0.042 
1.945 
0.911 


0.044 
2.089 
0.948 


0.034 
.068 
.502 


0.010 
2.021 
0.446 


-0.002 
-.144 
-.037 


0.043 
0. 143 
0.655 


0.042 
.072 
.613 


0.042 
0.084 
0.620 



-0.012 
-.007 


+ 0.001 
+ .071 
+.042 



Table 3 shows a comparison of total integrated 



absorption 



•/: 



yds, through the entire atmosphere 



calculated from actual radiosonde data and the total 
integrated absorption obtained from fits of the bi- 
exponential model (5 is the path length). There are 
three stations chosen, Balboa, Canal Zone (August), 
representing tropical conditions, Sopporo, Japan 
(August), representing temperate conditions, and 
Barrow, Alaska (February), representing arctic con- 
ditions. Each station is examined at three different 
frequencies: 6 Gc/s and 42 Gc/s on both sides of 
water vapor resonance, and the water vapor resonance 
frequency, 22.406 Gc/s. The initial elevation angle, 
0o, of the path is taken to be 90°, for purposes of 
simple integration of the exponential model. It is 
worthwhile to note that, even for the worst fits (wet 
exponential fit of Barrow, Alaska, which, as men- 
tioned before, fit best with a negative H w , thus account- 
ing for the negative values of Tw) the error in the total 
integrated absorption is not any larger than the other 
errors listed in table 3. 



3. Climatology of Surface Dry and Wet 
Absorptions 

A discussion of the equations for absorption in the 
6 to 45 Gc/s frequency range, and how these equations 
depend upon the commonly measured meteorological 
parameters of pressure, P, temperature, 7\ and water 
vapor density, p w , can be found in Bean [1964]. It 
is these formulas that were used for absorption com- 
putations on the profiles for the 49 stations mentioned 
earlier. The surface values of absorption from this 
49 station sample, were used as the basis for the con- 
tours drawn in figures 8, 9, 10, and 11. Figure 8 shows 
contours of ya, for a wavelength of 3 cm under winter- 
time (February) conditions, and figure 9 shows ya, 
contours for a 3 cm wavelength in the summertime 
(August). Regardless of their numerical value, for a 
particular contour map ya, and y w , should behave 
essentially the same at all frequencies in the range 
6 to 45 Gc/s, at least as far as meteorological and/or 
latitudinal variation is concerned. In figure 8, for 
example, a decided seasonal effect is noticed immedi- 
ately for this wintertime case when compared with the 
summertime picture of figure 9. The seasonal effect 
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is the fact that near the equator the values of yd, o 
remain essentially unchanged from summer to winter, 
but much higher values of yd, o are evident in arctic 
and upper northern hemisphere locations in the winter 
than in the summer, and vice versa for the antarctic 
and lower southern hemisphere locations. 

This variation is understandable if one recalls that 
yd, o is directly proportional to pressure, p , at the sur- 
face [Van Vleck, 1951]. Worldwide maps of average 
seasonal surface pressures show a decided analogy 
to figures 8 and 9 as regards the appearance of the 
contours. 

Figures 10 and 11 show maps of y w , o for wintertime 
(February) and summertime (August) conditions, 
where the values of y w , were obtained from the 49 
station sample discussed earlier. Because of the 
direct proportionality of y w> to surface water vapor 
density, p w , [Van Vleck, 1951], the most noticeable 
trend on figure 10 is the contenental-to-ocean change 
in y w , o especially in the northern hemisphere. This 
trend, no doubt, reflects the contrast of maritime air 
and the cold, dry air masses common to the interior 
of continents at this time of year. Surprisingly, at 
first, there seems not to be nearly so much of a similar 
effect in figure 11 in the antarctic region, but this is 
undoubtedly because there is so much more ocean 
between southern hemisphere continents than in the 
northern hemisphere. 

Figures 12 and 13 show the worldwide character of 
the dry absorption scale height, Hd, as determined 
from the 49 station sample. Again because of the 
independence of (3) with wavelength, climatic varia- 
tions on either map at one particular wavelength is 
representative of climatic variations for all wavelengths 
in the 6 to 45 Gc/s frequency range. Figure 12 repre- 
sents wintertime Hd (February) conditions for a wave- 
length of X = 1.339 cm, and figure 13 represents 
summertime H d (August) conditions at a wavelength 



of k — 3.0 cm. Note in general the latitudinal trend 
of the scale heights which closely follows the seasonal 
surface temperature trends for the world, again as 
expected from (3). 

Figures 14 and 15 represent the worldwide trend 
of the wet absorption scale height, H w , as determined 
from the 49 station sample. However, figures 2 and 
3 show that the essential important difference between 
Hd and H w in the 6 to 45 Gc/s frequency range is the 
maximum in the H w curve near water vapor resonance 
(22.4 Gc/s). This, of course, shows that H w is not 
independent of frequency. However, it can be shown 
that between 6 and 45 Gc/s: 



1 



— mfiipv) +/ 2 (P, To) +/ 3 (A, To) 



(6) 



where fi(p w ) is a function of water vapor density, p w , 
alone, and for that matter 



/.(p«.) 



Ho. 



(7) 



where H Pw is the scale height of an assumed exponen- 
tial fit of the form 



Pw = Po,w exp- 



UJ 



(8) 



to the water vapor density profile where po,w is the 
surface value of p w . The function /2(P, To) in (6) is 
a function of atmospheric pressure P, and surface 
temperature 7o, and / } (A, T ) is a function of wave- 
length and surface temperature. 
By differentiating (6), we obtain 



■7^A// u , = A/i + A/ 2 + A/3. 



(9) 




FIGURE 8. Characteristic worldwide wintertime (February) values FIGURE 9. Characteristic worldwide summertime (August) values 
of y d ,o. of y dj0 . 
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Figure 10. Characteristic worldwide wintertime (February) values Figure 12. Characteristic worldwide wintertime (February) values 
of y w , o- of H d . 




Figure 11. Characteristic worldwide summertime (August) values Figure 13. Characteristic worldwide summertime (August) values 
of y w , o. of H d . 



The A/ 3 term of (9) is unchanged, if temperature is 
unchanged, from location to location in figures 14 and 
15. Similarly, it would be unchanged on any other H w 
map from location to location when that map is plotted 
for a given frequency. For the normal atmospheric 
range of temperature conditions, A/ 3 generally contrib- 
utes less than 10 percent of (9). For the given condi- 
tions of figures 14 and 15, the variations of H w on any 
other map plotted for a different frequency should re- 
main about the same. Hence the seasonal climatic 
variations of H w for the frequency spectrum 6 to 45 
Gc/s are characterized in figures 14 and 15. It is 
interesting to note that for H Pw ^ 2.5 km, where H Pyo 
values are frequently clustered, f\ constitutes 85 per- 
cent or more of the total contribution to (6). 



Figure 14 indicates the wintertime (February) 
variation of H w (at k= 1.339 cm), and figure 15 indi- 
cates the summertime (August) variation of H w 
(at A. = 1.339 cm). The most important climatic feature 
to be noted on both maps is the continental-oceanic 
effect, but there are also effects depending upon the 
actual climatic character of a particular area. For 
example, in the summer there is a drop (fig. 15) in 
H w values towards the interior of North America. 
This is probably caused by the general continental 
tropical (dry) conditions that tend to prevail over the 
continent in the summer. The cross-hatched area in 
figure 14 represents regions of very poor fit of the biex- 
ponential model. These are the arctic cases, as indi- 
cated in table 4. Their presence is more severely felt 
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FIGURE 14. Characteristic worldwide wintertime (February) values 
o/H w . 

in the northern hemisphere in the winter (fig. 14) 
than in the summer (fig. 15) where, except at the poles, 
no appreciable arctic effect is present. 

Table 4. Ave rage correlation coefficients, 7, and standard errors 
of estimate S.E. of a fit of the form of yw^ywo exp(— z/H w ) at 
\= 1.339 cm to radiosonde data of several stations geographi- 
cally widely dispersed. 



Air mass 


Number of 

pieces of 

data averaged 


Mean surface^ 
temperature T {) 


Mean surface 

relative 

humidity RH 


- 


S.E. 


cA(7 , >-30°C) 

mT 

Combined cases: 

cP+mP + cT 


7 
55 

24 


•c 

- 15.9 
-1-21.8 

+ 6.1 


% 

80 
79 

76 


- 0.0065 
-0.9191 

-0.9098 


dBlkm 
±0.01796 
±0.03997 

±0.01942 



4. Climatology of Individual Biexponential 

Fits 

Since the oxygen absorption is easily fit with a 
descending exponential function of height with a 
correlation coefficient of —0.99 or better for all the 
cases tested in the 49-station sample, the most con- 
cern from location to location is the water vapor or 
wet absorption. One can easily find two different 
stations, even with quite similar geography and with 
quite similar surface meteorological conditions, yet 
exponential fits to their radiosonde profiles yield two 
very different values of scale heights, H w , correlation 
coefficients, r, and standard errors of estimate, S.E. 
Likewise, the same station under different synoptic 
weather conditions can have widely differing expo- 
nential fits. Unfortunately, unlike in the biexponential 
model of the refractive index [Bean, 1961] in the 6 to 
45 Gc/s region, the wet term of absorption usually 
tends to dominate the total integrated absorption (see 
figs. 8-11 for typical values), especially near water 
vapor resonance, X = 1.339 cm, thus making it relatively 



Figure 15. Characteristic worldwide summertime (August) values 
o/H w . 

imperative that the wet term exponential fit be closely 
examined climatologically. 

Perhaps the easiest climatological approach is 
classification of air masses according to source in 
standard meteorological fashion [Byers, 1944; Willett, 
1944] as: 

(a) Continental Arctic (cA) 

(b) Maritime Tropical (mT) 

(c) Maritime Polar (mP) 

(d) Continental Polar (cP) 

(e) Continental Tropical (cT). 

It would, of course, be necessary to examine the 
records of any particular location fully and perhaps, if 
frontal transitions are common, to consider mixed air 
masses [combinations of (a) through (e)] in order to 
determine the goodness of an exponential fit to the wet 
absorption. However, based on a study of 86 world- 
wide cases, table 4 helps to shed some light on the 
problem. 

The air masses were chosen from Willett's [1944] 
classification for the time of year considered. In 
table 4, one observes that the average correlation 
coefficient for cA cases is quite low. This is a uni- 
form tendency caused largely by the surface tempera- 
ture inversions noted in such air masses, and probably 
to some degree by the fact that radiosondes' humidity 
elements perform poorly at temperatures below °C 
[Wexler, 1949]. If the latter condition were corrected, 
the correlation coefficients in cA air masses might be 
considerably higher. It is also important to note 
that although the correlation coefficient is so low, the 
average standard error of estimate, S.E., is not very 
large. However, table 4 also shows that in spite of 
the much higher average correlation coefficient, an 
exponential model applied to the wet term for mT air 
mass cases might cause serious difficulties depending 
upon the application at hand. 
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5. Conclusions 

It has been seen that the fit of a biexponential 
model to atmospheric absorption with height is excep- 
tionally good for the dry term from 6 to 45 Gc/s and is 
generally good for the wet term in this frequency range 
except in cA cases where the wet term is small. An 
individual application requires the caution of the user 
and he should study his particular area climatologically 
and synoptically, depending upon the extent to which 
he plans application of the biexponential model. 



The authors wish to acknowledge the helpful advice, 
computer programs, and criticism of C. A. Samson, 
R. L. Abbott, E. R. Westwater, W. B. Sweezy, M. T. 
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ration of this paper. 
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